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Abstract 


The equilibrium configuration of a solid strange star in the final inspiral phase with another compact object is 
generally discussed, and the starquake-related issue is revisited, for a special purpose to understand the precursor 
emission of binary compact star merger events (e.g., that of GRB2112114A). As the binary system inspirals inward 
due to gravitational wave radiation, the ellipticity of the solid strangeon star increases due to the growing tidal field 
of its compact companion. Elastic energy is hence accumulated during the inspiral stage which might trigger a 
starquake before the merger when the energy exceeds a critical value. The energy released during such starquakes 
is calculated and compared to the precursor observation of GRB211211A. The result shows that the energy might 
be insufficient for binary strangeon-star case unless the entire solid strangeon star shatters, and hence favors a black 
hole-strangeon star scenario for GRB211211A. The timescale of the precursor as well as the frequency of the 
observed quasi-periodic-oscillation have also been discussed in the starquake model. 


Key words: stars: neutron — (stars:) gamma-ray burst: general 


1. Introduction 


It is well known that the puzzling nature of pulsar's interior is 
essentially relevant to the fundamental strong interaction at low- 
energy scale, the challenging non-perturbative quantum chromo- 
dynamics (or strong QCD Epelbaum et al. 2009), but this 
unknown state could be the first big problem to be solved in the era 
of multi-messenger astronomy (Baiotti 2019). Besides the conven- 
tional neutron star (NS) model, pulsars are proposed alternatively 
to be solid strange stars (Xu 2003) (or strangeon stars), and we are 
then developing a strangeon star model in order to understand 
extreme and mysterious events in astrophysics. Certainly, quakes 
can naturally occur on solid strangeon stars, and a giant quake 
model has already been proposed for the super-flares of isolated 
soft y-ray repeaters (Xu et al. 2006). Furthermore, tide-induced 
quakes in binary have also been discussed, appearing as so-called a 
sudden change in the tidal deformability at a certain breaking 
frequency of gravitational wave (GW) (Lai et al. 2019). Can tide- 
induced quakes be manifest in the electromagnetic (EM) wave? 
This is our focus here, and luckily, the precursor of GRB211211A 
could be a typical example. In the future, more such events, 
especially combined with the LVK-O4 observing run (e.g. 
Coupechoux et al. 2023), would surely be expected. 

Without doubt, the observation of GW170817 (The LIGO 
Scientific Collaboration & The Virgo Collaboration 2017) 


— black hole physics 


together with its EM counterparts GRB170817A and 
AT2017gfo (The LIGO Scientific Collaboration et al. 2017) has 
announced the birth of the multi-messenger astronomy era. This 
event has largely enriched our knowledge on the nature of short 
gamma-ray bursts (sGRB) (Narayan et al. 1992; LIGO Scientific 
Collaboration et al. 2017), the state of matter at supranuclear 
densities (Bauswein et al. 2017; Margalit & Metzger 2017; 
Rezzolla et al. 2018; Ruiz et al. 2018; Kiuchi et al. 2019; Shibata 
et al. 2019) as well as the origin of heavy elements in the 
Universe (Eichler et al. 1989; Abbott et al. 2017). The EM 
counterparts have been detected in almost every band from radio 
to gamma-ray, however, these observations all happen during the 
post-merger phase. As an implementation, EM signals prior to 
the merger (i.e., the precursor observation Troja et al. 2010), if 
detected, could significantly improve our understanding of the 
properties of the merging objects, as well as improve the 
detection and localization of the following GW signal. 
Interestingly, 2021 December 11th, a very peculiar gamma- 
ray burst (GRB) has been detected by Fermi/GBM (Mangan 
et al. 2021), Swift/BAT(D’Ai et al. 2021) and Insight- 
HXMT/HE (GRB211211A) (Zhang et al. 2021). An excess in 
optical/near-infrared has been identified, the multi-band 
properties of which are quite similar to those of AT2017gfo 
(Rastinejad et al. 2022; Xiao et al. 2022). Together with the 
non-detection of a supernova at the GRB location, this GRB is 
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suggested to be associated with merger event involving an NS, 
though the duration of the main emission is relatively longer 
(~8 s) compared with typical sGRBs. More intriguingly, a fast 
rising and exponentially decaying precursor has been observed 
approximately 1s prior to the main emission. The precursor 
lasts for ~0.2s and a quasi-periodic-oscillation (QPO) with 
frequency ~22 Hz has been identified in it (Troja et al. 2022; 
Xiao et al. 2022). 

Following this interesting observation, various models have 
been suggested to explain the precursor of GRB211211A. 
Previous force-free simulations of NS magnetosphere have 
shown that the interaction of the magnetic fields of two 
inspiralling-in NSs could produce Poynting flux strong enough 
to be observed as precursor emissions (Carrasco & Shibata 
2020). If a magnetar is involved in the merger event, a 
catastrophic flare of the magnetar during the inspiral phase 
could also be the source of the precursor (Xiao et al. 2022; 
Zhang et al. 2022). In addition, a resonant shattering of the 
solid crust of the merging NSs is also invoked to explain the 
observation, with certain demands on the NS spin and magnetic 
field (Suvorov et al. 2022). 

In this paper, we come up with a starquake model to explain 
the precursor of GRB211211A based on a solid strangeon star 
scenario, in which the energy budget could be satisfied 
regardless of the spin and magnetic field strength of the 
merging NSs. In this model, the equilibrium configuration of 
the solid compact star changes as the binary gets closer in the 
inspiral stage and tidal field of the companion becomes 
stronger. Stress will be accumulated as the elastic structure 
resists the change in the configuration. Eventually, the elastic 
strain might exceed a critical value before merger and the solid 
structure of the star cracks (i.e., a starquake happens). Thus, a 
precursor will be triggered by energy released during the 
starquake and the reconfiguration of the star. 

The paper will be organized as the follows: in Section 2, we 
will introduce the configuration of a solid strange star in the 
tidal field of its companion; a quantitative comparison with the 
observation will be made in Section 3; future observational 
prospects of this scenario will be discussed in Section 4. 


2. The Model 


2.1. Equilibrium Configuration of a Solid Strange Star in 
Close Binary 


The equilibrium configuration of a solid strange star in the 
tidal field of its companion is determined by the bulk energy of 
the star E;ota, Which consists of several parts, 


Egg = Ex + E, + E, + Eana, (1) 


in which E; is the kinetic energy of the star when rotation is 
considered, E, stands for the additional gravitational binding 
energy of the star, Ej, is the elastic energy accumulated and E, 
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is the energy possessed by the star due to the tidal field of its 
companion. 

It is believed that an NS spins slowly before merger due to 
magnetic dipole radiation during the long inspiral stage. 
Moreover, it has been shown that the tidal interaction is 
insufficient to synchronize the NS spin before merger (Bildsten 
& Cutler 1992). Hence, it is reasonable to assume that both E; 
and the change in E, are negligible before merger. 

E, and E, are related to gravity, and it is necessary to define a 
zero energy configuration as a reference. In the following, we 
will assume the spherically symmetric configuration to possess 
energy E>. In binary systems, however, the shape of the star 
will no longer be spherical and it is useful to define the shape of 
the star by the parameter of reduced ellipticity 


NIS 
lo 


€ 


(2) 
in which J is the moment of inertia of the star with respect to 
the axis of the deformation and Jp is that of the spherical star 
with the same baryonic mass. It is worth noting that, in the case 
of a rotating star, the star is flattened with respect to the rotation 
axis and e will be positive. For tide-induced deformation, the 
star will be prolonged along the direction where the companion 
is. Hence, e will be negative. For the consideration of this 
calculation, rotation effect is negligible as mentioned before 
and the principal axis is chosen to be the direction of the 
companion. For incompressible star and small ellipticity, € is 
related to the geometrical eccentricity of the star (e) as 
e--e /3. With this definition, one can obtain 


M. 3 
E, By = Eo + Age? + 1A.) e. (3) 


in which M and R are the mass and radius of the star, M, is the 
mass of its companion, D is the separation of the binary system, 
A, and A, are coefficients related to gravitational binding 
energy which depend on the density distribution of the star. In 
the case of an incompressible star (which is a good 
approximation for strange stars), the following relation holds 


3 GM? 
ee 4 
s 25 R e 
and 
2 
PNE (5) 


in which & is a dimensionless constant of order of unity. Ideally 
K = 3/5, but may vary due to the fact that the tide-induced 
deformation may not be perfectly along the direction of the 
companion or density structure of the star. For simplicity, we 
will use the notation A — GM? /R in the expressions below and 
A ^ 1 x 1054 erg for typical NS value M—1.4M. and 
R=10 km. 
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The elastic energy Ej, is related to difference between the 
reference ellipticity e; when the star solidified (for example, 
when the star suffers the previous starquake) and the current 
value e. According to Hooke's Law, one can obtain 


Eala = Ble = eo» (6) 


in which B — SHV, p and V are the shear modulus and volume 
of the solid strange star, respectively. 
Collecting all the ingredients we have 


3 M.( RY 
Evota = E + sae? + nAtte( ee Be- a 2 (7 
total ot as im ( o' TM 


and the ellipticity of equilibrium configuration could be 
obtained by minimizing the total energy which requires 
OE ota /O€ = 0. Thus, we could obtain the reduced ellipticity 


25A u R ) 50B 


K—-——_ — | — — €. 
6A + 50B M\D 6A + 50B 


For a purely fluid star (or for a solid star when starquake 
happens, to reconfigure itself), the equilibrium ellipticity is then 


25 M.( RJ 
Ce => -h . 9 
" SA) (9) 


It is worth noting that the solid star tends to resist from being 
deformed due to the accumulation of elastic energy. Hence, if 
the star was at its fluid equilibrium configuration before the 
previous solidification, then at any time during the later 
evolution, we should have € e,44 (note that € is negative). 
Additionally, the tidal interaction becomes stronger as the 
binary gets closer during the inspiral stage, the ellipticity of the 
star would increase with time, so we could obtain the following 
inequality which holds during the entire inspiral stage 


E0 > € > €eq,fl- (10) 


(8) 


Previous calculation and observations indicate that the shear 
modulus of a strange star j lies in the range between 10?? and 
10°% erg cm? (Xu 2003; Zhou et al. 2004). Within this range, 
the value of B would be much smaller than A. Combining the 
inequality Equation (10), one could obtain that the second term 
in Equation (8) is negligible compared with the first term. 
Therefore, we will omit the second term in Equation (8) in the 
calculations below and assume the equilibrium configuration of 
a solid strange star satisfies 


25A M,(RY 
€eq,so — —K . (11) 
6A + 508 M \ D 


2.2. The Starquake Model 


Equation (11) indicates how the ellipticity of the solid 
strange star increases as the separation of the binary shrinks 
during the inspiral stage, due to the dissipation of angular 
momentum and energy through GW radiation. Elastic energy 
gradually increases as the shape of the solid star changes. 
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Depending on the microscopical model of the solid strange star, 
the maximum stress of the solid structure might be reached 
before the merger happens and hence a starquake takes place. 
During the starquake, the ellipticity of star tends to migrate 
from its solid equilibrium configuration (€eq,so) to its fluid case 
(€eq,a). Therefore, not only the elastic energy but also the 
change in E, + E, will be released during the starquake as they 
are related to the ellipticity. This starquake scenario is quite 
similar to the starquake model of pulsar glitches (Baym & 
Pines 1971; Zhou et al. 2004, 2014; Lai et al. 2018). In the 
latter case, the accumulation of the elastic energy (i.e., the 
change in the ellipticity) is due to the change of pulsar spin, 
which results from the magnetic dipole radiation of pulsars. 

We will first calculate the change in E, + E, before and after 
the starquake. The change of the reduced ellipticity is the 
difference between €,4,5 and Eeq,so: 


3 
"CS AA). (12) 
6A + 50B 6 M\D 


According to the feasible range of u, B is at least two orders of 
magnitude smaller than A and hence we have de « 0.1e,g. 
Therefore, we could estimate the change in E,-- E, by an 
expansion with respect to e: 


O(E, + Ej) 
Oe 


However, e,4 is obtained by requiring O(E, + E,)/de=0. 
Consequently, the non-vanishing parts of E, + E, are terms of 
6€ and higher order ones and hence are negligible. Similar 
conclusion is found for the case of starquake scenario of pulsar 
glitches (Zhou et al. 2014). As a result, the major contribution 
of energy release during such starquake is the release of the 
accumulated elastic energy before the starquake: 


Eg, = B(€eq,so is eo)? 
25A M. R ) 
= B| -k 
6A + 50B M \ D 
25 M,( R Y] 

gos (14) 

6 M\Do 
in which Do is the separation of the binary system when the 
solid strange star experienced its previous solidification. As the 
precursor happened very close to the merger (i.e., when the 
orbital separation shrinks rapidly), it is reasonable to make the 


following assumption Do > D. In addition, we remind that 
A> B, then the elastic energy is approximately 


3 2 
px ae [ (5) | | (15) 


óc = €eq,fl — €eq,so = — 


(E, + Ej) Lunes (13) 


It is easy to verify that the total change in E, + E, is indeed 
~10% of Eaa at most and we will focus on Equation (15) when 
comparing with the observation in the next section. 
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3. Results 
3.1. Energy Budget 


The origin of GRB211211A is still under debate. Scenarios 
including binary neutron star (BNS) merger (Kunert et al. 
2023), black hole-neutron star (BH-NS) merger (Zhu et al. 
2022) as well as neutron star-white dwarf merger (Yang et al. 
2022) have been proposed. In this section, we will focus on 
BNS and BH-NS merger scenarios and test their feasibility and 
parameter space in producing the observed precursor within the 
starquake scenario. 

The major difference between BNS and BH-NS merger 
scenarios is their mass ratio M./M and the orbital separation D 
when the starquake happens. For the BNS scenario, the most 
mass-asymmetric BNS system observed and the mass of which 
is precisely measured has mass ratio q = 0.78 (Ferdman et al. 
2020). GW190425 was measured to have an even smaller mass 
ratio of 0.7, whereas it is still uncertain that whether its heavier 
component is an NS or a BH (The LIGO Scientific Collaboration 
et al. 2020). Even though there are equation of state (EoS) 
models for solid strange stars which could reach maximum mass 
as high as above 3 Mo (Lai & Xu 2009), this could merely push 
the possible range of q to about 0.5 (or, M./M — 2.0 when we 
consider the lighter star suffers the starquake). 

For the BH-NS merger case, the range of the binary mass 
ratio could be much wider. GW observations indicate that 
Mgu/Mws could be as small as 2.0, if heavier components of 
those mass gap events are indeed BHs. The upper limit for 
Magu/Mws for the event GRB211211A should be constrained 
by the observation of associated kilonova. If the mass of the 
BH is too much larger than NS, the ejected mass from the NS 
would be insufficient for powering a kilonova even with an 
extreme BH spin. According to previous studies, we put 10.0 as 
a maximum possible value for Mgy/Mnys (Kyutoku et al. 2015; 
Kawaguchi et al. 2016). In the analysis below, we will assume 
M./M to be in the range of [2.0, 10.0] for the BH-NS merger 
scenario. 

The orbital separation D when the starquake happens could 
be implied by the time when the precursor happens. For 
GRB211211A, the time between the precursor and the main 
burst is approximately one second, which means the starquake 
happens less than one second before merger happens as there 
might be time delay between the actual merger and the time 
when jet is launched (for instance, for the BNS merger 
GW170817, the time delay between the merger and sGRB is 
approximately 1.7 s). For BNS system this indicates an orbital 
separation D < ~ 100 km according to Figure 1. For the BH- 
NS system the possible range of D is larger (i.e., could be as 
large as ~300 km), due not only to a wider range of the mass 
range, but also to BH spin which could affect the dynamics of 
the final inspiral phase. Particularly, extreme spin is necessary 
for significant mass ejection to take place for BH-NS system 
with large mass ratio. 
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Figure 1. The evolution of the orbital separation D prior to the merger (the 
merger time is calibrated to t= 0 in this figure) for BNS and BH-NS cases, 
obtained by effective-one-body (EOB) calculations (Damour & Nagar 2014; 
Nagar et al. 2016, 2018, 2020a, 2020b; Riemenschneider et al. 2021). The 
green solid and dashed curves on the top represent the BH-NS merger scenario 
with M./M = 10 and BH spin of —0.8 (i.e, anti-aligned with the orbital 
angular momentum) and 0.8, respectively. The black solid and dashed curves 
on the bottom stand for the BNS merger cases with M,/M = 1.0 and 2.0, 
respectively. 


For the BNS case, inserting the value of 50 km for D and 
unity for «x, one can obtain the total elastic energy contained in 
the entire star as the following 


Ea, = 2.3 x 109 erg] — —— 
( 10?^ erg cm? 


2 -6 
x (=x) (sz ) ; (16) 
1.0 50 km 
For the observation of GRB211211A, the energy released in 
the precursor in the form of electromagnetic emission is 
~7.7 x 105 erg which means the total energy released should 
be larger than this value. Therefore, the energy budget is very 
tense for explaining the precursor in the BNS case, unless the 
solid structure of the entire star shatters and the starquake has to 
happen at a closer range with the jet launching being delayed 
after the merger. Alternatively, one can count on a very mass- 
asymmetric merger with M./M as large as 2.0. Nevertheless, 
even in this case, the total elastic energy is not enough if D is 
larger than 100 km when the starquake happens. 
The detailed result is shown in Figure 2. Even if we allow for 
a global starquake during which the elastic energy of the entire 
star is released and all converted into precursor EM emissions, 
the orbital separation at the moment of this starquake has to be 
smaller than ^60 km for an equal mass BNS merger. This 
upper limit on D could be relaxed to ^75 km if we consider 
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Figure 2. Constraints on the binary parameters according to the precursor 
observation of GRB211211A in the BNS case. Three dashed lines labeled by 
100%, 50% and 10% represent the combination of the mass ratio M./M and 
orbital separation D when the starquake happens with which the corresponding 
percentage of the total elastic energy has to be released in order to explain the 
energy of the precursor observation. The gray shaded region is excluded since 
more than 100% of the total energy budget is required (and hence insufficient). 
If we consider that partial failure of the solid structure (i.e., less than 10% of the 
elastic energy is released) is more realistic, the possible parameter space is very 
narrow unless the starquake happens extremely close to the merger (D smaller 
than 50 km). 


extreme mass ratio cases. Nevertheless, as mentioned above, 
such global starquake is less likely to happen and a partial 
starquake which releases 10% of the elastic energy contained in 
the star is not possible unless the starquake happens at a binary 
separation less than 50 km and with large mass ratio. 
Compared with the BNS case, the BH-NS merger scenario is 
more favored considering the energy budget due to a larger 
possible M./M. Ten percents of the total elastic energy 
converting into EM emission would be enough for the observed 
precursor luminosity. According to previous researches, it is 
shown that elastic energy is not uniformly accumulated 
throughout the solid star when deformation is induced. Partial 
failure in the solid structure is also demonstrated to be more 
reasonable to explain the observation of pulsar glitches (Lai 
et al. 2018). In our case, the deformation is induced by the tidal 
field of the companion of a solid strange star. Consequently, the 
stress is expected to be the largest near the surface of the star, 
especially near the polar and equatorial region. Thus a small 
fraction of the solid star shatters and releases the elastic energy 
in this region is a more natural scenario. Indeed, there would be 
larger uncertainties in inferring the orbital separation by the 
time before merger for the BH-NS case. Therefore, we treat D 
as a free parameter and explore the possible combination for D 
and M../M for different amount of elastic energy required. 
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Figure 3. Constraints on the binary parameters according to the precursor 
observation of GRB211211A in the BH-NS case. Three dashed lines are 
plotted in the same way as in Figure 2. The darker bottom left gray shaded 
region is excluded according to the energy budget whereas the lighter gray 
shaded region on the top right is inside the ISCO of the BH companion 
(assuming that the mass of the NS is 1 Mo and the BH is a Schwarzschild BH) 
and hence is excluded. Nevertheless, the ISCO of a spinning BH can be much 
smaller than that of a Schwarzschild BH and parameters in this region could as 
well be possible, as extreme spin is essential for sufficient mass ejection in the 
case of large mass ratio BH-NS mergers. 


The results for the allowed parameter space of the BH-NS 
case are shown in Figure 3. As the larger allowed mass ratio 
benefits the accumulation of elastic energy, the separation of 
the binary could be as large as ~100 km when the starquake 
happens if we allow for a 100% release and conversion of the 
elastic energy. There is still plenty of possible parameter space 
even if we consider a partial (10%) elastic energy release for 
orbital separation D ~ 60 km. For even smaller fraction (such 
as 196 case), the binary separation at which the starquake 
happens lies inside the innermost stable circular orbit (ISCO) of 
a Schwarzschild BH with the same mass of the BH companion. 
However, as mentioned above, large spin is needed for large 
mass ratio (M/M ~ 5) BH-NS merger to produce a kilonova 
and the ISCO radius of a Kerr BH could be much smaller than 
that of the Schwarzschild case. Consequently, it is still possible 
for a tiny starquake to happen at such small separation and 
account for the precursor observation of GRB211211A. 


3.2. Sequential Starquakes During the Insipiral 


The above discussion is based on the assumption that the 
deformation induced is homogeneous in the entire star and only 
one starquake happens prior to the merger, which could either 
be a global quake or a local one. Nevertheless, the later in the 
inspiral stage, the faster the binary orbit shrinks (see Figure 1) 
and hence the elastic energy is accumulated much more 
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Figure 4. An illustration of the sequential starquake model as discussed in 
Section 3.2. The first starquake happens at fo with binary orbital separation Do. 
After that, sequential starquakes take place with increasing depth inside the 
solid strange star, as the binary gets closer. At a later time (for instance, 
to + 0.2 s), only the matter inside a sphere with radius Ry could maintain its 
solid structure, with Rp approximately satisfies the following relation: 
Rj/ Dy — R/Do. The gray shaded region illustrates the part of star of for which 
the elastic energy could be released during this sequential starquake. 


rapidly. Moreover, the ellipticity induced by the tidal field 
increases as (R/Dy which suggests the deformation of the star 
is larger in its outer part for a realistic star (1.e., which could not 
be incompressible) at a fixed separation. With these considera- 
tions, it is more likely that several sequential starquakes take 
place during the last of the inspiral phase, if the first starquake 
is a partial one. However, in this case elastic energy is partially 
released and we focus on the BH-NS scenario in this section, 
due to its larger parameter space. 

Our sequential starquake model is demonstrated in Figure 4. 
The first starquake takes place at a separation of Do for the 
binary at time tọ. As the deformation is largest in the outer part, 
only the surface part of the star suffers the quake (for instance, 
the region from R — 6R to R, in which ôR is the depth of this 
starquake). At a later time t; the binary would be separated at a 
smaller distance of D; and a shell of the star at a radius of Ry 
would experience the same deformation as those at the surface 
of the star at the separation of Do, in which Rp is simply 
determined by the geometrical relation R;/D;= R/Dyo. Thus, 
the matter at Ry meets its limit for a failure in its elastic 
structure. In this case, several individual starquakes take place 
during the time fo to t; from the surface of the star to the depth 

The observations of the precursor of GRB211211A indicate 
that the energy is released in several individual bursts with 
decaying peak amplitude in a time span of ~0.2 s, which is also 
the origin of the 22 Hz QPO. Such a feature could be well 
explained by a sequential starquake model: those later 
starquakes happen deeper inside the star and consequently 
the energy release tends to be diffused more slowly and hence 
the peak flux observed would be lower. 
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Figure 5. Constraint on the binary separation when the last quake takes place 
in the sequential quakes scenario. The mass ratio and BH spin is M./M = 10 
and s — 0.8 for this model. With the green dashed curve we plotted the fraction 
of the volume of the star which suffers starquakes in the duration of 0.2 s 
a — (Dr /Do)), which is equal to the fraction of elastic energy released, 
assuming a uniform shear modulus. The black curve shows the ratio between 
the observed precursor energy and the total elastic energy (Eops/Eeta). The 
intersection point then sets up an upper limit for the parameter Dy. 


In such a model, the fraction of elastic energy released 
during the sequential starquakes over a certain duration 
depends on the binary separation of the first/last quake. As a 
result, in the sequential starquakes scenario, there will be 
additional constraint on the parameter space. The fraction of the 
elastic energy released in 0.2 s could be obtained as 


(17) 


3 
Evelease DN Vauake —1- (2) 
Eala Votal Do 


in which Eea also depends on Dy as in Equation (16). 
Therefore, given the mass ratio and BH spin, Do could be 
determined by the EOB model and hence Ejejease could be 
obtained. The fact that the energy released has to be greater 
than the observed value then sets a constraint on Dj. 

The result is shown in Figures 5 and 6, in which we have 
chosen the case of M,/M — 10 and 3, respectively. As 
mentioned above, for larger mass ratio case, large BH spin is 
also essential (Kyutoku et al. 2015; Kawaguchi et al. 2016) and 
we have chosen s= 0.8 for the dimensionless spin for the 
M./M = 10 case. For smaller mass ratio, we have explored 
three different BH spin (s = 0, 0.4 and 0.8) parameters to verify 
its impact on our results. Our result shows that a constraint is 
not very sensitive on the BH spin and upper limit for D; is 
approximately 120km for M./M=10 and 75km for 
M./M —3, which are all much larger than the ISCO radius. 
It is worth noting that, if one assumes a certain energy 
conversion efficiency (7) from the released elastic energy to the 
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Figure 6. The same as Figure 5 but for a mass ratio of M,/M = 3. In this case 
we have also explored different values of BH spin to check its influence on our 
analysis, the green dashed curve, red dashed curve and blue dashed curve stand 
for the cases of s — 0, 0.4 and 0.8, respectively. 


EM emission, a tighter constraint could be made by requiring 


n= Eobs/Eela > (18) 
— (DgJDo) 
which could be read from the result figures once a value of 7) is 
assumed. 
In this model, the rising timescale of the precursor could also 
be understood. The change in the ellipticity of the star after the 
first quake happens approximately in a free falling timescale as 


27R 


tise ^9 —. 
rise GM./ Do 


The timescale is larger for smaller mass ratio case as M./Do 
becomes smaller according to the analysis above. In the 
M. M = 3.0 case, the above formula yields a rising timescale 
of ~0.8 ms, which is consistent with the observations. 

In addition, the QPO frequency in the observation of the 
precursor of GRB211211A could also be naturally understood 
as the number of quakes during the entire process within 0.2 s. 
However, the number of quakes and the energy released in each 


(19) 


quake in the entire sequential process depend on the thickness 
of the shells which suffer the failure of the solid structure 
during each quake. The thickness of the shattered shells could 
be very sensitive to properties which are difficult to predict 
such as crystal defects in the solid structure which might be 
caused during previous quakes. Hence, the observational 
properties (i.e., the QPO frequencies) could be totally different 
even for sources with similar binary parameters. 
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3.3. Quake-induced Oscillation? 


The enormous energy released during a tide-induced giant 
quake may set the entire star into vibration, producing 
oscillations at frequencies determined by the structure and 
elastic properties of the solid strange star. Typically, non-radial 
oscillations of solid stars include spheroidal and toroidal 
modes. Toroidal (torsional) modes are a type of oscillation that 
maintains the star’s shape. They are purely shear oscillations. 
QPOs following giant flares in soft gamma-ray repeaters and 
anomalous X-ray pulsars indicate the shear motion of the star 
after a giant quake (Duncan 1998; Watts & Strohmayer 2007). 
The spheroidal oscillations refer to waves that change the shape 
of the star, involving both radial and tangential displacements. 
These can be easily excited in a tidally induced starquake. Both 
types of oscillations have been observed in large earthquakes 
(Benioff et al. 1961; Park et al. 2005). Furthermore, these 
modes couple more easily with external magnetic fields than 
modes coming from the deep interior of the star. For these 
reasons, we expect that the tide-induced torsional modes may 
couple with Alfvén waves along the magnetospheric field lines, 
resulting in the QPO in the precursor of GRB211211A. 

We denote (£, £5, Eg) as the displacement amplitudes in 
spherical coordinate. To estimate the frequencies of toroidal 
and spheroidal modes, we model the solid strange star as a 
homogeneous, isotropic elastic sphere with uniform density 
and shear modulus. 

The toroidal oscillations are divergence free with no radial 
components. For a particular eigenmode denoted by / and m, 
the separation of variable for the displacements takes the 
form (McDermott et al. 1988) 

m Wr) OY, iwt 


OY,, i 
se ; Q——W — 1 e't 20 
Šo sinÜü OQ i Šo e) 00 E eo 


where Y,,,(0, @) are spherical harmonics. Inserting Equation (20) 
into the shear wave equation, we obtain the following 
equation for the radial eigenfunction W(r) (McDermott et al. 
1988) 


Id + 1) 
H r? 


dr? r dr 


2 2 
dw 2dw E LE QD 


where p and u are the star's density and shear modulus 
respectively. The torsional modes are referred by the notation 
nt} Where n is the overtone number of radial nodes in the 
eigenfunction W(r). In this paper, we focus on the nodeless 
vibrations with n = 0. By applying the boundary condition of 
vanishing surface horizontal traction, we can solve the 
eigenvalue problem analytically (Lamb 1881), obtaining the 
eigenfrequency 


f(ot) = Cr (22) 
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Here v = (u/ p)" ? is the shear speed, C; is a constant depending 


on l. The fundamental mode of torsional oscillation is of; 
because ofp has zero displacement and of, describes a constant 
azimuthal twist of the entire star, which does not exist for free 
oscillation. 

We denote the spheroidal mode as os;. The separation of 
variable for the displacement can be written as 


Yin V(r) OY 
es = 


E E o 
Er = U(r)Y, So = VG) 50" sind Op (23) 


Combined with perturbation in the gravitational potential, we 
can obtain the systematic differential equations governing the 
oscillation. We refer the readers to McDermott et al. (1988), 
Crossley (1975), Alterman et al. (1959) for those equations. To 
simplify the calculations and represent the eigenfreqeuency in a 
simple relation to the shear speed, we use the long wavelength 
approximation. The eigenfrequency of nodeless modes is 


[204 + DE- DI? v 
2m R 


The oscillation of order /— 1 does not exists, and the 
fundamental mode of the global nodeless spheroidal oscillation 
is 052. 

Both the strong and the electromagnetic interactions are 
involved in characterizing the shear modulus, j4 of solid 
strange matter. One may have p ~ 10?^ erg cm ^? if the strong 
force dominates the shear mode Xu (2003). In the catastrophic 
process of the tide-induced giant quake, the star may be 
fractured as a whole and release enormous energy. The 
modulus could be decreased after the fracture. 

From Equations (22) and (24), we obtain the fundamental 
frequencies for the torsional and spheroidal modes 


f (os) = (24) 


1 


» u i 
x 
= H i 
f (oS2) = 61.7 ncs . (26) 


Here we take M = 1.4 Mo and R = 10 km. There is one thing to 
be clarified, here we are discussing about the quake-induced 
oscillation modes, which are modes propagating in a solid star 
after the starquake. Therefore, the shear modulus involved in 
the calculation has to be smaller than that before the starquake 
which is applied in the calculations in Section 3. One can 
notice that the mode frequency is on the same order of the QPO 
frequency observed in GRB211211A if the shear modulus of 
the solid strange star decreases to the order of 
1099 10? ere cm ?. In the mean time, the required shear 
modulus after the starquake is much smaller (i.e., 3 orders of 
magnitude smaller) than those applied before the starquake, and 
hence has only negligible influence on the calculation of the 
energy budget. As a conclusion, the model is capable of 
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Figure 7. The relation between eigenfrequencies of the fundamental modes of 
the torsional and the spheroidal oscillations. We show the cases with / — 2, 
3, 4. 


explaining the energy budget and QPO frequencies at the same 
time. In Figure 7, we show the fundamental frequencies with 
1 — 2, 3, 4 as functions of the shear modulus. 

In reality, the frequency of a mode will be shifted due to the 
effects of gravitational redshift and the Doppler effect. The 
gravitational redshift factor caused by the NS can be expressed 
as 


2GM Y? 
s = (1 - G ) 1, (27) 


where M and R represent the mass and radius of the NS, 
respectively. For example, assuming an NS with a mass of 
1.4 Mo and a radius of 10 km, the gravitational redshift factor 
is 0.306. 

To estimate the redshift factor caused by the binary system, 
we neglect the spin of the companion and crudely treat the NS 
as a test particle. We assume that the observer is located in the 
orbital plane and at an orbital phase of ¢ = 0 when the NS is 
closest to the observer. We select several special points to 
discuss the frequency shift. At o = 0, 7, the redshift factor is 
given by: 


3GM. y? 
Do ) 1, (28) 
where M, represents the total mass of the binary system, and D 
represents the separation of the binary system. This equation 
takes into account both the gravitational redshift and the 
transverse Doppler effect caused by the motion of the NS in the 
binary system. 
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At $ = 7/2, the redshift factor for the binary system is 
given by 


For a binary system with M, = 10 Mc, D = 100 km, M = 1.4 Mo, 
and R=1l0km, the redshift — factor  Zbinary(@ = 0) = 
Zpinary(P x 1/2) — 0.34, Zpinary(P mi 1/2) =0.95, and Zbinary(® a 
—7/2) — —0.27. The redshift Zvinary is largest at o = 7/2. The 
frequency is blueshifted at o = —7/2 since the Doppler effect 
dominates over the gravitational redshift. 

The relation between the observed frequency fobs and the 
mode frequency in the local frame of the NS, f, can be 
approximated by: 


Tas = (1 + Zn)! + Zbinary) Tf. (30) 
The largest redshift corresponds to fos, (o = 7/2) = 0.393f. 


4. Conclusions and Discussions 


In this paper, we introduced the scenario of a tidal induced 
starquake for a solid strange star prior to the merger with a 
compact companion. As the binary inspirals closer, the tidal field 
will gradually deform the solid strange star, resulting in the 
accumulation of the elastic energy. A giant starquake might be 
triggered when the binary is close enough and the stress exceeds a 
critical value. This scenario is quite similar to the starquake model 
of pulsar glitches, in which the elastic energy is accumulated as 
the star deforms due to the spinning down process. 

In particular, we demonstrated that such giant quake before 
merger could release sufficient energy to explain the precursor 
observation of GRB211211A. In the sequential quakes model, 
which we consider only part of the star suffers the starquake, 
the released energy is still enough at a separation larger than the 
ISCO radius. Moreover, the torsional mode frequencies are 
estimated by regarding the entire solid strange star as a 
homogeneous, isotropic elastic spheroid with uniform density 
and shear modulus. The result depends on the final shear 
modulus when starquake happens (i.e., after the failure of the 
solid structure takes place) and could be consistent with the 
observation with a reasonable range of the final shear modulus, 
without contradicting the energy budget calculations. 

Although we have obtained sufficient energy budget to 
account for the total energy release and oscillation frequencies 
consistent with the QPO frequency of the observation, the 
mechanism for converting those energy into EM radiation 
(particularly in the hard X-ray and gamma-ray band) and the 


8 In fact, since we do not assume any details about the EoS model or mass- 


radius relation of the compact star, any model that could provide the required 
shear modulus would be able to explain the observations within our scenario. 
However, it is unlikely that NSs which is believed to have a superfluid core 
could accumulate to the required elastic energy budget before a starquake takes 
place. 
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way to connect asteroseismology frequency with QPO is more 
complicated and is often not discussed in details in previous 
studies of GRB211211A (e.g., Suvorov et al. 2022) and QPOs of 
soft gamma-ray repeaters (Duncan 1998). Nevertheless, there are 
simulations indicating that accelerated charged particles in the 
magnetosphere of an orbiting NS might be capable of producing 
emissions in the X-ray and even gamma-ray bands (Carrasco & 
Shibata 2020; Carrasco et al 2021). Meanwhile, stellar 
oscillation is shown to be related to pulsar activities by affecting 
the magnetosphere (Lin et al. 2015). Future investigations will 
be needed to connect these ingredients together and to better 
understand the entire scenario. 

Our result favors a BH-NS scenario compared with the BNS 
case, due to much larger parameter space for the consideration of 
the energy budget. The BH-NS scenario is also more consistent 
with the event rates. SGRBs with confirmed precursor observa- 
tion are quite rare, only as BH-NS binary systems are much 
fewer compared with BNS systems. We will estimate the event 
rates in the future and try to verify this argument. 

We are expecting to test the model presented in this paper in 
the future, especially by combining the observations of gravita- 
tional and electro-magnetic signals. The LVK-O4 observing run 
(e.g., Coupechoux et al. 2023) would start, and China's mega- 
facilities, especially the Gravitational wave high-energy Electro- 
magnetic Counterpart All-sky Monitor (GECAM) (Li et al. 2020), 
the Hard X-ray Modulation Telescope (HXMT) (Li et al. 2018), 
as well as the planned Einstein Probe (EP) (Yuan et al. 2022), 
would work. The model we proposed might then be soon 
falsified, but, conversely, would show its particular reasoning 
style in the coming years. 
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